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Abstract 

The  interaction  of  intense,  femtosecond  laser  pulses  with  a  dielectric 
medium  is  examined  using  a  numerical  simulation.  The  simulation  uses 
the  ID  electromagnetic  wave  equation  to  model  laser  pulse  propagation. 
In  addition,  it  includes  multi-photon  ionization,  electron  attachment, 
ohmic  heating  of  free  electrons,  and  temperature  dependent  collisional 
ionization.  Laser  pulses  considered  in  this  study  are  characterized  by  peak 
intensities  ~  1012  -  1014  W/cm2  and  pulse  durations  ~  10-100  fsec.  These 
laser  pulses  interacting  with  fused  silica  are  shown  to  produce  above- 
critical  plasma  densities  and  electron  energy  densities  sufficient  to  attain 
experimentally  measured  damage  thresholds.  Significant  transmission  of 
laser  energy  is  observed  even  in  cases  where  the  peak  plasma  density  is 
above  the  critical  density  for  reflection.  A  damage  fluence  based  on 
absorbed  laser  energy  is  calculated  for  various  pulse  durations.  The 
calculated  damage  fluence  threshold  is  found  to  be  consistent  with  recent 
experimental  results. 
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I.  Introduction 

The  interaction  of  high-intensity  femtosecond  laser  pulses  with  dielectric 
materials  involves  a  number  of  physical  mechanisms  that  occur  on  time  scales 
comparable  with  both  the  laser  frequency  and  pulse  duration,  and  spatial  scales  that  can 
be  shorter  than  the  laser  wavelength.  For  example,  the  electron  collision  time  can  be 
comparable  with  the  optical  laser  period  and  the  width  of  the  plasma  layer  generated  at 
the  dielectric  surface  can  be  much  less  than  an  optical  wavelength  [1].  In  addition,  the 
collisional  ionization  time  of  the  dielectric  surface  is  comparable  with  the  pulse  duration. 
In  these  situations  an  accurate  description  of  laser  pulse  propagation  is  necessary  to 
describe  the  interaction  and  subsequent  transmission  of  laser  energy  into  the  dielectric. 

A  high-intensity  (>1013  W/cm2),  ultrashort  (-100  fsec)  laser  pulse  impinging  on  a 
dielectric  surface  initially  generates  free  electrons  (plasma)  through  multiphoton  or 
tunneling  ionization.  For  the  femtosecond  time  scales  involved,  these  electrons  gain 
energy  from  the  laser  field  through  collisional  heating.  As  the  electron  temperature 
increases  the  collisional  ionization  rate  also  increases.  The  plasma  density  generated  on 
the  dielectric  surface  can  be  greater  than  the  critical  density  for  reflecting  the  laser  pulse. 
The  photoionization  and  collisional  ionization  process  for  femtosecond  laser  pulses  has 
been  studied  extensively  both  theoretically  [1-5]  and  experimentally  [1,6-10]. 

Theoretical  descriptions  of  the  collisional  ionization  process  range  in  complexity  and 
include  fully  kinetic  models  [4],  models  which  make  the  flux-doubling  assumption  [1, 

10],  and  multiple  rate  equation  descriptions  [5]. 

In  these  previous  studies,  laser  propagation  is  treated  in  an  approximate  fashion, 
i.e.,  the  electric  field  in  the  dielectric  is  taken  to  be  evanescent  and  the  ionization 
equations  are  solved  purely  in  the  time  domain  at  a  given  spatial  position  for  a  given 
temporal  laser  profile.  Stuart,  et  al.,  [1]  considered  a  limited  laser  propagation  model 
involving  the  steady  state  form  of  the  transmission  coefficient  to  calculate  the  transmitted 
laser  pulse  shape  and  energy  deposition  depth.  That  study,  using  1053  nm  and  526  nm 
wavelength,  sub-picosecond  pulses,  predicted  above  critical  plasma  densities  and  plasma 
layer  widths  smaller  than  the  laser  wavelength.  For  a  laser  pulse  with  1053  nm 
wavelength,  -200  fsec  duration,  and  3.5  J/cm2  fluence,  i.e.,  above  the  damage  threshold, 
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Stuart’s  propagation  model  predicted  significant  absorption  of  laser  energy  (~60%) 
within  1  pm  of  the  surface. 

Here,  we  utilize  a  fully  time-dependent  theoretical  model  to  simulate  the 
interaction  of  a  high-intensty,  femtosecond  laser  pulse  with  a  dielectric.  We  model  the 
propagation  of  the  laser  in  ID  using  the  electromagnetic  wave  equation  coupled  to  an 
equation  for  the  induced  plasma  current.  The  plasma  is  generated  by  multiphoton  and 
collisional  (avalanche)  ionization.  The  collisional  ionization  process  is  described  by  a 
reduced  kinetic  model  which  allows  us  to  write  a  temperature  dependent  ionization  rate 
in  closed  form.  The  plasma  electrons  are  heated  by  the  laser  field  via  collisions  with  ions 
and  neutrals.  For  femtosecond  laser  pulses,  the  absorbed  electromagnetic  energy  of  the 
laser  is  mainly  converted  to  kinetic  energy  of  free  electrons.  The  hot  plasma  layer  formed 
on  the  surface  of  the  dielectric  is  typically  much  smaller  in  width  than  an  optical 
wavelength.  The  heated  electrons  transfer  their  energy  to  the  lattice  on  a  time  scale  much 
longer  than  the  laser  pulse  duration  [1]. 

The  model  equations  are  solved  numerically  to  determine  laser  pulse 
transmission,  laser  energy  absorption,  damage  thresholds,  and  plasma  density  and 
temperature  within  the  dielectric.  Our  results  show  much  greater  transmission  compared 
with  a  purely  steady  state  calculation.  The  enhanced  transmission  is  due  to  the  fact  that 
the  critical  density  for  reflection  is  generated  only  after  a  significant  part  of  the  pulse  has 
crossed  the  surface  boundary.  To  a  lesser  degree,  the  high  collisionality  of  the  plasma 
layer  also  enhances  transmission.  The  transmission  of  laser  energy  through  the  plasma 
layer  has  applications  to  the  operation  and  design  of  plasma-based  switches  and  mirrors 
[]  and  possibly  electronic  countermeasures. 

In  this  study  we  consider  laser  intensities  of  1013  -  10 14  W/cm2,  and  pulse 
durations  from  tens  to  hundreds  of  femtoseconds,  i.e.,  pulses  that  can  be  propagated 
many  tens  of  meters  in  air  [12],  We  limit  our  study  to  parameter  regimes  for  which 
thermal  and  particle  diffusion,  nonlinear  self-focusing,  and  group  velocity  dispersion  can 
neglected.  We  note  that  others  have  used  full  PIC  simulations  to  model  laser  propagation 
in  dielectrics,  although  the  parameter  regimes  investigated  are  different  from  what  is 
considered  here  [13].  Experimentally,  the  transmission  of  focused,  intense  femtosecond 
laser  pulses  through  fused  silica  was  observed  by  Von  der  Linde  and  Schueler  [14]. 
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Section  II  of  the  paper  discusses  the  theoretical  model.  Section  III  presents  the 
results  and  compares  with  experimental  data  from  single  pulse  interactions  with  fused 
silica.  Conclusions  are  presented  in  Sec.  IV. 


II.  Theoretical  Model 

We  consider  the  situation  in  which  a  laser  pulse  propagating  along  the  z-axis  is 
incident  on  a  dielectric  which  is  located  at  z>0,  as  shown  in  Fig.  1 .  The  laser  intensity  is 
high  enough  that  it  can  singly  ionize  the  dielectric  and  heat  the  resulting  plasma.  The 
laser  electric  field,  E  =  E(z,t)x,  and  magnetic  field,  B  =  B(z,t)  y,  are  determined  by  the 
Maxwell  equations 

r)  F  r)  /? 

n20(z, 0—  =  -c~  -  An[Jx(z,t)  +  ,  (la) 

ot  oz 


dt  dz 


(lb) 


where  the  source  current  Jx  denotes  the  linear  plasma  response  and  Jmpi  is  an  effective 

current  which  accounts  for  laser  energy  depletion  due  to  multi-photon  ionization  [15]. 

The  refractive  index  is  given  by  a  low  frequency  Lorentz  model,  i.e., 

n\  (z,t)  =  1  +  a>l(z,t)/  Q] ,  where  the  quantity  Q0  is  a  constant  that  determines  the 

magnitude  of  the  refractive  index  in  the  absence  of  ionization, 

a>n(z,t )  =  -^4  nq2nn(z,t)l  m  ,  and  nn(z,t )  is  the  density  of  neutrals.  In  writing  the 


refractive  index  this  way,  dispersion  is  neglected. 

The  linear  transverse  plasma  current  is  given  by 


+ 


a2P(z>t) 

An 


E, 


(2) 


where  cop(z,t)  =  -yj 4nq2ne ( z,t )/ m  is  the  electron  plasma  frequency,  ne  is  the  electron 

density.  Equation  (2)  is  exact  for  the  case  of  a  cold  plasma  in  which  free  electrons  are 
created  with  zero  mean  momentum  [15].  Viscosity  effects  have  also  been  neglected.  The 
electron  collision  frequency  ve  =  vei  +  ven  is  the  sum  of  the  electron-neutral  and 
electron-ion  collision  frequencies  which  are  given  by 
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V„  [sec-1]  =  2xl0-7  n,  [cm"3]  T,m  [eV] ,  (3a) 

v„[sec-|]  =  2.91xl0-6Z«,[cm-3]A(»„7;)r,-3,![eV];  (3b) 

respectively,  where  A (ne,Te)  is  the  Coulomb  logarithm  and  Te  is  the  electron 
temperature  [16].  For  this  study,  we  limit  Ato  values  greater  than  2. 

The  plasma  generation  is  described  by 
d  yi 

~J7  =  Vn>  nn  +  Vcol  (Te  T,rap  ,  (4) 

where  vm  and  vcol  are  the  multi-photon  and  collisional  ionization  rates,  respectively.  For 
this  study,  we  use  the  6-photon  ionization  rate  for  the  interaction  of  0.8  pm  radiation  with 
fused  silica  (ionization  energy  Wjon  =  9  eV ), 

vjsec’1]  (5) 

where  cr6  =2xl013  cm-3  psec-1  (cm2/TW)6  [17]  is  the  Keldysh  6-photon  ionization 
coefficient,  nn0  is  the  neutral  density  in  the  absence  of  ionization,  and  I  is  the  slowly 

varying,  time  averaged  laser  intensity.  It  is  assumed  that  only  single  ionization  of 
neutrals  occurs.  This  assumption  will  not  significantly  affect  the  transmission  of  the  laser 
pulse  since  all  transmission  occurs  before  the  saturation  of  ionization,  i.e.,  the  critical 
density  for  reflection  is  approximately  three  times  lower  than  the  maximum  ionization 
density. 

The  exponential  decay  of  electron  density  given  by  the  third  term  on  the  right  side 
of  Eq.  (4),  has  been  observed  experimentally  with  a  characteristic  decay  time, 
ztrap  =  150  fsec  [18].  Physically,  the  overall  decay  in  electron  density  is  due  to  a  number 
of  processes,  e.g.,  recombination,  attachment,  diffusion. 

The  current  density  Jmpi  can  be  derived  by  equating  the  time  rate  of  change  of 

laser  energy  density  due  to  multiphoton  ionization  with  the  product  JmpiE .  The  result  is 
[15] 

JmPi  =  Wwnvm(E)nJE,  (6) 

where  Wi0„  is  the  ionization  energy. 
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jto 

The  collisional  ionization  rate  can  be  written  as  vcol  =  v0(W)  fe(W )  dW  , 

•0 


where  fe  is  the  electron  energy  distribution  function,  and  v0  (W)  =  a0(W  IWion  -  l)2  for 
W  >  Wton  and  zero  for  W  <  Wjon  [2].  The  rate  constant  a0  is  a  particular  to  the  material 
under  consideration  and  can  be  empirically  determined.  For  fused  silica,  a0  =1.5  fsec'1 
[1,19].  We  assume  a  Maxwellian  distribution  for  the  electron  velocity,  i.e., 

F(v)  =  (nmvth)'3  exp(-v2  !v]h) .  The  corresponding  energy  distribution  function  is 


f.m  =  2 


/  \>/2 
‘  W  ' 


exp  (-W/WJ.. 


where  Wlh  =  mv\  12  =  2>kBTe  12  and  fe(W)  is  normalized  such  that  J°  fe(W)  dW  =  1 . 
Using  this  distribution  function,  the  collisional  ionization  rate  is  given  by 


Vco,(/?)  = 


2«o  nn 


n 


no 


15 


>3/2 


2 


exp(-l/^)  +  2^^15 


P2  -3(3  +  \  Erfc(l  /  4P ) 

4  ) 


(7) 

* 

where  /?  =  Wth  /Wion  and  Erfc(x)  =  1  -  (2  / 4ft)  Jexp(-t2 ) dl  denotes  the  complementary 

o 


error  function. 

For  the  parameter  regime  considered  here,  thermal  conduction  can  be  neglected 
and  the  electron  thermal  energy  satisfies  the  equation 


w 

dt 


JXE  -  T  Wionvco,  ne 


(8) 


The  first  term  on  the  right  side  of  Eq.  (8)  represents  ohmic  heating  of  the  electrons 
while  the  second  term  represents  the  energy  lost  by  free  electrons  due  to 
collisional  ionization.  The  constant  T  takes  into  account  the  fact  that  due  to  other 
electronic  excitations,  it  usually  requires  energies  greater  than  Wion  to  ionize  by 


collisions.  For  example,  in  gas  discharges,  2  <  T  <  5  for  electron  energies 
comparable  with  the  ionization  energy  [20].  Here,  we  take  Y  =  2  as  a  canonical 
value.  Eq.  (8)  implicitly  assumes  that  electrons  generated  by  multiphoton 
ionization  are  created  with  zero  kinetic  energy,  i.e.,  above  threshold  ionization  is 
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neglected.  It  is  also  assumed  that  in  the  electron  decay  process,  an  electron  that  is 
captured  transfers  its  kinetic  energy  to  another  electron.  This  is  apparent  from  the 
lack  of  a  cooling  term  in  Eq.  (8)  that  is  proportional  to  the  electron  decay  rate 

l/Tr rap- 

III.  Results 

Equations  (1),  (2),  (4),  and  (8)  are  solved  numerically  using  a  finite  difference, 
staggered  leap-frog  approach.  In  what  follows  we  consider  laser  pulses  with  wavelength 
A  =  0.8  pm ,  durations  ranging  from  30-400  fsec,  and  fluences  <  6  J/cm2.  The  temporal 
profile  of  the  laser  electric  field  with  carrier  frequency  co0  at  a  fixed  z-position  in  vacuum 
(  z  <  0 )  is  given  by 

E(z  <  0,  t)  =  E0  sin[;r t/rL\ cos[n>0 ( t  -  ( rL  / 2))]  (9) 

for  0  <  t  <  zL  and  zero  otherwise.  The  full-width-half-max  of  intensity  is  given  by 

TFWHU  =  /  2  . 

The  laser  pulse  is  made  to  interact  with  a  dielectric  located  in  the  region  z  >  0 . 
The  dielectric  is  taken  to  be  fused  silica,  i.e.,  refractive  index  n0  =1.5,  ionization  energy 
Wmn  =9  eV ,  and  neutral  density  nn0  =  2  x  10 22  cm-3 .  The  ionization  rate  constant  is 
aQ  =1.5  fsec-1  [1,19].  For  the  assumed  wavelength  and  ionization  energy,  multiphoton 

ionization  is  a  6  photon  process  where  the  ionization  rate  is  described  by  Eq.  (5).  Figure 
2  plots  the  multi-photon  ionization  rate  as  a  function  of  laser  intensity  (Fig.  2a)  and  the 
collisional  ionization  rate  as  a  function  of  J3  =  Wth  /  Wion  (Fig.  2b)  for  the  assumed 
dielectric  parameters. 

Steady  State  Solution 

Laser  pulses  that  are  sufficiently  intense  will  ionize  the  dielectric  surface, 
producing  a  thin  plasma  layer  that  can  reflect  and  absorb  laser  energy.  Before  we  solve 
the  fully  time-dependent  problem  including  ionization,  we  first  consider  the  steady  state 
interaction  of  a  laser  with  a  pre-existing  thin,  highly  collisional  plasma  layer  on  the 
surface  of  a  dielectric  and  calculate  the  transmission  of  the  laser  into  the  dielectric.  We 
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take  the  region  z  <  0  to  be  vacuum,  the  region  from  0  <  z  <  L  to  contain  a  pre-ionized 
dielectric  with  a  given  electron  temperature,  and  the  region  z  >  L  to  contain  a  dielectric  of 
refractive  index  n0 .  Ionization  and  heating  by  the  laser  pulse  is  neglected  in  this 
simplified  steady  state  model.  It  is  assumed  that  the  regions  z  <  0  and  0  <  z  <  L  contain 
right  and  left  going  waves,  while  only  a  transmitted  (right-going)  wave  exists  in  the 
region  z  >  L.  The  electron  collision  frequency  is  calculated  according  to  Eqs.  (3).  If  the 
plasma  density  and  temperature  within  the  region  0  <  z  <  L  are  uniform,  the  transmitted 
intensity  for  z  >  L  can  be  calculated  analytically.  The  curve  in  Fig.  3  shows  the 
transmitted  intensity  (normalized  to  the  incident  intensity)  versus  the  ratio  cop  /  ©0 .  In 

generating  this  curve,  we  use  parameters  characteristic  of  the  fully  time  dependent 
simulations  presented  later,  i.e,  electron  temperature  of  10  eV,  n0  =  1.5 ,  X  =  0.8  pm, 

and  L  =  1/10.  It  is  seen  that  for  maximum  ionization,  i.e.,  cop  =  3.4  <y0 ,  the  steady  state 

calculation  predicts  a  fractional  transmission  of-  3  x  10“2 . 

In  the  following  section,  we  show  that  significantly  higher  transmission  is  observed  when 
the  fully  time-dependent  equations  are  solved.  To  benchmark  the  time-dependent 
numerical  code,  we  simulate  the  steady  state  interaction  described  above.  The  results  of 
the  simulation  are  denoted  by  the  dots  in  Fig.  3.  Excellent  agreement  is  obtained  with  the 
analytic  result. 

Time  Dependent  Simulations 

We  now  consider  the  fully  time  dependent  numerical  solution  of  Eqs.  (1),  (2),  (4), 
and  (8).  Figure  4  shows  a  simulation  in  which  a  laser  pulse  with  rL  =100  fsec ,  and 
fluence  3.2  J/cm2  (/0  =  l.lxl 014  W/cm2)  impinges  on  fused  silica.  Figure  4a  shows  the 

spatial  profile  of  the  laser  electric  field  before  the  interaction  with  the  dielectric  located  in 
the  region  z  >  0.  After  interaction  with  the  dielectric  surface,  reflected  and  transmitted 
pulses  are  observed,  as  shown  in  Fig.  4b.  A  plasma  layer  is  generated  with  a  maximum 
plasma  frequency  cop  =  3.4  <a0 ,  which  corresponds  to  full  ionization  of  the  dielectric 

surface.  The  half-maximum  width  of  the  plasma  layer  is  -  0.2  pm.  A  portion  of  the 
plasma  layer  near  the  surface  is  strongly  heated  by  the  laser  pulse.  For  this  particular 
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simulation  the  electron  temperature  reaches  a  peak  value  of  ~  10  eV.  The  half-maximum 
width  of  the  heated  electron  layer  is  ~  0.05  pm. 

Even  when  the  plasma  density  is  above  critical,  i.e.,  cop  >  co0 ,  approximately  30% 

of  the  incident  laser  energy  is  transmitted  through  the  surface  as  a  propagating  pulse.  The 
transmitted  pulse  is  steepened  at  the  leading  edge  and  elongated  at  the  back.  It  continues 
to  ionize  the  dielectric  as  it  propagates.  However,  its  intensity  is  sufficiently  low  that  the 
plasma  density  remains  well  below  the  critical  value.  The  spatial  profile  of  the  reflected 
pulse  shows  a  gradual  rise  in  amplitude  at  the  front  followed  by  a  sharper  rise  towards  the 
middle  of  the  pulse,  indicating  a  sudden  increase  in  the  reflectivity  of  the  surface. 

Figure  5  plots  the  normalized  plasma  density  at  the  surface  versus  time  for  the 
simulation  shown  in  Fig.  4.  The  laser  pulse  is  located  in  the  region  0  <  t  <  100  fsec . 

The  plasma  density  reaches  the  critical  value  ~  40  fsec  into  the  interaction.  Hence,  the 
sharp  rise  in  the  reflected  pulse  amplitude  seen  in  Fig.  4  is  coincident  with  the  plasma 
density  reaching  the  critical  value.  For  comparison,  Fig.  5  also  plots  the  plasma  density 
versus  time  for  purely  multi-photon  ionization  ( a0  =  0 ).  Comparison  of  the  two  curves 

shows  that  collisional  ionization  increases  the  peak  plasma  density  by  30%.  Hence, 
while  photoionization  is  the  dominant  ionization  mechanism,  collisional  ionization  also 
plays  an  important  role. 

Figure  6  shows  the  peak  plasma  density  (6a),  fractional  transmitted  fluence  (6b), 
and  fractional  reflected  fluence  (6c)  from  several  simulations  utilizing  different  pulse 
durations  and  intensities.  Transmitted  fluence  is  measured  at  position  z  =  1  pm.  At 
smaller  incident  fluences,  for  which  ionization  is  negligible,  i.e.,  cop  «  co0 ,  the  pulse  is 

reflected  and  transmitted  according  to  Snell’s  Law.  At  higher  fluences,  when  ionization 
becomes  important,  more  of  the  laser  energy  is  reflected  and  the  transmission  decreases. 
However,  greater  than  70%  transmission  of  fluence  is  observed  even  when  the  plasma 
density  reaches  the  critical  density.  Greater  than  40%  transmission  can  still  occur  when 
the  dielectric  is  fully  ionized  (cop  =  3.4 o)0). 

Strong  absorption  of  laser  energy  occurs  near  the  surface.  The  amount  of  laser 
energy  absorbed  depends  on  the  duration  and  intensity  of  the  incident  laser  pulse.  In  Fig. 
7,  we  plot  the  fluence  absorbed  within  1  pm  of  the  surface  versus  incident  fluence  for 
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various  pulse  lengths.  The  absorbed  fluence  is  calculated  as  the  incident  fluence  minus 
the  sum  of  the  reflected  and  transmitted  fluences,  where  the  transmitted  fluence  is 
measured  at  z  =  lpm.  The  parameters  are  the  same  as  in  Fig.  6.  It  is  seen  that  for  smaller 
incident  fluences,  cop  «  co0 ,  and  there  is  no  significant  absorption  of  laser  energy.  As  the 

laser  intensity  is  increased  for  a  fixed  pulse  duration,  ionization  occurs  and  the  dielectric 
surface  becomes  highly  absorptive  when  cop  >  co0 .  Figure  6  shows  that  for  incident 

fluence  <  2  J/cm2,  short  pulses  deposit  energy  into  the  dielectric  more  efficiently  than 
long  pulses. 

Since  the  transition  from  transparency  to  opacity  occurs  very  rapidly,  we  have 
also  tried  slightly  different  pulse  shapes  to  see  if  the  reflection  and  transmission  would  be 
significantly  affected.  We  have  run  simulations  with  smoother  sine-squared  field  profiles 
and  found  only  small  changes  to  the  transmission  and  reflection  coefficients.  For 
example,  for  the  case  of  a  100  fsec  pulse  the  reflection  coefficient  increases  by  ~10% 
with  the  smoother  pulse  shape. 

Damage  Threshold 

Experimentally,  damage  has  been  characterized  as  an  observable  modification  of 
the  surface  when  viewed  through  a  Nomarski  microscope  [1,8].  Experimental  results  also 
suggest  that  a  minimum  energy  density  must  be  deposited  within  the  material  in  order  to 
cause  observable  damage.  The  energy  density  barrier  for  fused  silica  is  Ub  =  54  kJ/cm3 
[10]. 

Using  our  simulation,  we  can  define  a  damage  threshold  based  directly  on  energy 
absorption.  Let  fc  (z0 )  denote  the  total  (cumulative)  laser  fluence  (J/cm2)  absorbed 

within  a  layer  0  <  z  <  z0 .  The  energy  density  absorbed  by  the  dielectric  is  then  given  by 
Uabs  =  /c(z0)/z0.  We  define  the  damage  fluence,  fd ,  as  the  minimum  incident  laser 
fluence  required  to  achieve  Uabs  =Ub.  Figure  8  plots  /c(z)  and  Uabs(z )  for  the  same 
parameters  as  Fig.  4.  The  dependence  of  fabs  on  z  indicates  that  most  of  the  laser 

absorption  takes  place  within  ~0.5  pm  of  the  surface.  The  maximum  absorbed  energy 
density  is  -180  kJ/cm3,  which  is  well  above  what  is  required  for  damage. 
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Other  analyses  have  defined  the  damage  fluence  as  the  minimum  incident  fluence 
required  to  create  a  plasma  with  a  density  such  that  cop  =  a >0 ,  arguing  that  the  plasma 

becomes  highly  absorptive  when  the  critical  density  is  attained.  However,  because  the 
plasma  density  can  reach  the  critical  density  after  a  significant  portion  of  the  laser  pulse 
has  already  passed  through  the  plasma  layer,  the  condition  cop  =  coQ  may  be  insufficient 

for  femtosecond  laser  pulses.  For  example,  a  100  fsec  laser  pulse  with  an  incident 
fluence  of  ~1.3  J/cm  can  create  a  plasma  with  the  critical  density  (see  Fig.  6a). 

However,  the  maximum  absorbed  energy  density  for  this  interaction  is  ~  2.5  kJ/cm3, 
which  is  much  less  than  Ub . 

Figure  9  plots  the  damage  fluence,  fd ,  calculated  using  two  different  definitions, 

versus  pulse  duration.  The  solid  dots  denote  the  damage  fluence  calculated  according  to 
our  definition.  The  open  circles  denote  the  fluence  required  to  create  a  plasma  with 
cop  =  co0,  i.e.,  the  damage  fluence  as  defined  in  Refs.  [1,6-8].  The  damage  fluence  is 
seen  to  increase  with  pulse  duration  from  -1.4  J/cm2  to  4.4  J/cm2  over  the  range 
50  fsec<  tL  <  400  fsec .  The  damage  fluence  denoted  by  the  open  circles  is  lower  by 
approximately  40%  over  this  range  of  pulse  durations.  For  example,  for  a  pulse  duration 
of  tl  -  100  fsec ,  it  requires  -2.2  J/cm  of  incident  laser  fluence  to  cause  54  kJ/cm  of 
laser  energy  to  be  absorbed  within  the  dielectric. 

IV.  Conclusions 

The  interaction  of  intense  femtosecond  pulses  with  a  dielectric  surface  has  been 
numerically  modeled  in  ID  using  the  electromagnetic  wave  equation,  rate  equations  for 
the  electron  density,  and  an  equation  describing  electron  energy  density.  The  rate 
equations  contain  multiphoton  ionization,  collisional  ionization,  and  electron  trapping. 
The  equation  for  the  electron  energy  density  describes  ohmic  heating  and  cooling  by 
collisional  ionization.  Because  the  simulation  solves  the  time-dependent  electromagnetic 
wave  equation,  the  transmission,  reflection,  and  absorption  of  the  laser  pulse  by  a  thin, 
highly  collisional  plasma  layer  is  modeled  for  the  first  time  in  this  parameter  regime  in  a 
fully  self-consistent  manner.  Significant  differences  are  observed  between  the  fully  time- 
dependent  simulation  and  the  approximate  steady  state  model.  It  is  seen  that,  due  to  the 
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time  dependence  of  the  ionization  process  a  significant  transmission  of  laser  energy  is 
possible  even  when  the  plasma  density  is  above  the  critical  value.  This  is  because  the 
density  can  exceeds  the  critical  value  towards  the  trailing  end  of  the  laser  pulse.  For 
example,  consider  the  transmission  of  a  pulse  with  tL  =100  fsec  and  fluence  of  3.2 
J/cm2.  The  interaction  of  this  pulse  with  the  dielectric  fully  ionizes  the  surface  and  heats 
the  plasma  to  ~10  eV.  From  Fig.  6b  the  fractional  transmission  of  this  pulse  is  ~  25% 
according  to  the  time-dependent  solution  while  the  corresponding  steady  state  calculation 
(Fig.  3)  predicts  a  fractional  transmission  of  only  3%. 

The  simulations  also  show  that,  for  parameters  characteristic  of  recent 
experiments  using  fused  silica,  the  absorbed  energy  density  is  sufficient  to  cause 
observable  damage.  Our  damage  fluences  are  in  good  agreement  with  the  experiments  of 
Refs.  [7,8].  We  note,  however,  that  the  damage  fluence  in  those  experiments  were 
defined  as  the  incident  fluence  required  to  cause  observable  discoloration  of  the  surface 
when  observed  through  a  Nomarski  microscope.  For  our  simulations,  the  damage 
threshold  is  defined  as  the  incident  fluence  required  to  deposit  an  energy  density  of  54 
kJ/cm3  within  the  material  [10].  This  energy  density  is  sufficient  to  break  a  significant 
number  of  the  molecular  bonds  near  the  surface. 
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Figure  1 :  Schematic  illustration  of  the  numerical  model  geometry. 
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d)p/co0 

Figure  3:  Steady  state  transmitted  intensity  normalized  to  the  incident 
intensity  versus  normalized  plasma  frequency  for  X  =  0.8  pm  radiation 

impinging  on  a  dielectric  ( n0  =1.5)  with  a  preformed  surface  plasma 
layer  of  width  L-  X!  10  and  electron  temperature  of  10  eV. 

Ionization  and  heating  by  the  laser  pulse  are  neglected.  Solid  curve 
denotes  analytic  result,  dots  denote  simulation  results. 
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Figure  4:  (a)  Initial  laser  field  E  vs.  z  before  the 
interaction,  t  =  0  defines  the  time  at  which  the  leading 
edge  of  the  pulse  reaches  the  dielectric  surface  located  at 
z  =  0.  (b)  Electric  field  of  the  transmitted  and  reflected 
pulses  t  =  122  fsec  after  the  interaction.  Darker  curve 
denotes  cop  / a>0  vs.  z.  Pulse  duration  is  rL  =100  fsec, 
with  fluence  3.2  J/cm2. 
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Figure  5:  Peak  value  of  cop  /  o)0  vs.  time  for  cases  of 
multiphoton  ionization  only  (  a0  =  0  ),  and  both  multiphoton 
and  collisional  ionization  ( a0  =150  fsec“').  Pulse  duration 

r\ 

is  tl  =100  fsec ,  with  fluence  3.2  J/cm  . 
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Figure  6:  (a)  Normalized  plasma  frequency,  (b)  transmitted 
fluence,  frl  fn  and  (c)  reflected  fluence  fref  /  fi  versus  incident 

fluence,  fj  for  pulse  durations  zL  =  50, 100,  200,  and  400  fsec . 
Transmitted  fluence  was  recorded  at  position  z  =  1  pm. 
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Figure  7:  Fluence  absorbed  within  the  dielectric  from  z  =  0  to  z  =  1  pm  for  the  same 
parameters  as  Fig.  6. 
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Figure  8:  Cummulative  absorbed  laser  fluence,  fc ,  and 
absorbed  energy  density,  Uabs  versus  z  at  time  t  =  122  fsec 
for  the  same  parameters  as  Fig.  4. 
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Figure  9:  Damage  fluence,  fd ,  versus  pulse  duration,  rL  for 

fused  silica.  The  damage  fluence  denoted  by  the  solid  dots  is 
defined  as  the  minimum  incident  fluence  required  to  absorb  54 
kJ/cm3  of  laser  energy  within  the  dielectric.  The  damage  fluence 
denoted  by  open  circles  is  defined  as  the  incident  fluence 
required  to  create  a  plasma  with  cop  =  a>0. 
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